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Abstract 



We show how the isospin vector and axialvector current spectral functions 
pv,3 and PA.3 can be used to determine in leading chiral order the low energy 

(3/2) (3/2) 

constants and . This is accomplished by matching the Oper- 

ator Product Expansion to the dispersive analysis of vacuum polarization 
functions. The data for the evaluation of these dispersive integrals has been 
recently enhanced by the ALEPH measurement of spectral functions in tau de- 
cay, and we update our previous phenomenological determination. Our calcu- 
lation yields in the NDR renormalization scheme and at renormalization scale 
fi = 2 GeV the values ''^^ = 0.55±0.07±0.10 and = 1.11±0.16±0.23 

for the quark mass values m<j + m = 0.1 GeV. 
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I. INTRODUCTION 



The recent KTeV and NA48 findings that e'/e ~ 20 ■ 10 ^ raise the important question 
whether a value so large can be consistent with Standard Model expectations. [|[] One of the 
key quantities upon which the Standard Model prediction is based is the B-factor -Bg . 

/o /OA 

In this paper, we work in the chiral limit to obtain an analytic expression for Bg (and 
also for B) ). Our results take the form of sum rules involving the difference pv — Pa of 
isospin-one spectral functions. 

The constants B^^^'^^ {p) and B^^"^^ (p) are defined in terms of the matrix elements 

(27r|Qf/^)|K), = Bf'/'\p) (2vr|Qf/^)|ir);- = 7,8) , (1) 

where p is the renormalization scale, Q'^^'^^ and Q^^'^^ are the A/ = 3/2 electroweak penguin 
operators 

= SaT^dh {ubT'^Ua - dkT%) + SaT^UhUtT% , (2) 

a, b are color labels and Fl = 7^(1+75), Fj^ = 7^(1—75). Analogous B-factors are defined for 
the other weak operators. The most important contributions to e'/ e are the matrix elements 
of the penguin operator Qg and the electroweak penguin operator Qg . As expressed in 
terms of the B-factors this is seen in the approximate numerical relation. 



e' 



10 X 10" 



, (3) 



when evaluated aX p = 2 GeV in the MS'-NDR renormalization scheme. Alternatively, in 
terms of the operator matrix elements themselves one has 



- = 10 X 10"^ 
e 



-3.1 GeV-3 • (Q6)o - 0.51 GeV-^ ■ {Qt'^^)2] (4) 
where 

(Q6)o = ((7r7r),=o|Q6|i^°) and {Qf"\ = ((7r7r),=2| , (5) 

again defined in the NDR scheme at /i = 2 GeV. 

A kaon-to-pion weak matrix element can be analyzed by chiral methods and expressed 
as an expansion in momenta and quark masses. Here we shall calculate the leading 
term in such an expansion, valid in the limit of exact chiral symmetry where the Goldstone 
bosons become massless. In fact, our analysis yields values for the relevant K-io-n matrix 
elements themselves. However, we shall transcribe this information into the equivalent form 
of B-factors in order to express our results in the more conventional form and thus allow 
comparison with other techniques. When considered in the chiral limit, the content of Eq. (|l]) 
reduces to 
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where Oi, are the local four-quark operators 
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^1 = (lli^lT^ ^T^'tt^ 



T-3 
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?7AtA — g g7''A — g - g7M75A 



(7) 



In the above, q = u,d, s, is a Pauli (flavor) matrix, {A^} are the Gell Mann color matrices 
and the subscripts on Oi, refer to the color carried by their currents.^ A chiral evaluation 
of Bj^^'^^ and _Bg is thus seen (c/Eq. (|^)) to depend upon {Oi)^^ and {O^)^. The tttt 
matrix elements are likewise recoverable in the chiral limit from the vacuum matrix elements. 
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The rest of this paper describes a calculational procedure to obtain analytic expressions for 
the vacuum matrix elements. 

In Section II, we show how to extract {Oi) ^ and (Os)^ from the isospin-one vector and 
axialvector correlators by deriving sum rules for {Oi)^ and {Os)^i in a momentum cutoff 
scheme. This renormalization is especially well suited for comparing theory to experiment. 
In Section III we introduce a four-quark operator Oas=i, distinct from the familiar nonlep- 
tonic hamiltonian TCAs=iy whose kaon-to-pion matrix element is related to (Ci)^ and (Cs)/^ 
in the chiral limit. We demonstrate consistency of this information with the sum rules of 
Section II. Section IV describes a procedure for obtaining (Oi)^ and (Cs)/^ in MS renormal- 
ization, which is commonly used in lattice-theoretic simulations. Our final numerical results 
and concluding statements appear respectively in Sections V and VI. 



II. DISPERSIVE ANALYSIS OF VACUUM POLARIZATION FUNCTIONS 

In seeking values for {Oi)^ and (Cs)/^, it is natural to study the vacuum polarization 
functions as these are also defined in terms of vacuum matrix elements of four-quark opera- 
tors (but generally not all at the same spacetime point). Thus we consider the combination 
11^,3 — n^^3 (the subscript '3' denotes the isospin flavor), 

t J d'x e'^-{0\T{Vr{x)V,%0) - A^,{x)Am) |0) 

= (gV-gV'^)(ny,3-n^,3)(g')-gVn2k?') • (9) 



Throughout we denote vacuum expectation values as (0|O|0)^ = (O)/^. 
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Associated with this correlator is the difference of spectral functions pv^ — Pa,3, 



[ny,3 - n^,3] (Q' 



1 



ds 



[PV,3 - PA,3] (s) 



(10) 



where = — g^. In writing this spectral relation, we have made use of the first and second 
Weinberg sum rules 0, which are both valid in the chiral limit 

Due to the complexity of QCD, there exist no analytic expressions for the correlators 
and spectral functions that are valid over the entire energy domain. However, some crucial 
information is available. At low energies, pv,3 and pA,3 are determined from r-lepton decays 
and from e~^e~ scattering. As one proceeds from the resonance region of nonperturbative 
physics to larger energies, the effect of individual channels becomes indistinguishable and 
perturbative QCD (pQCD) becomes operative. The boundary between nonperturbative and 
perturbative regions defines a scale A ~ 2 — 3 GeV. In the pQCD domain, the leading-log 
behavior of (Jlv,3 ~ ^a,3){Q^) is given by 



Q\Uv,3 - ^A,3){Q') ~ 27r{asOs), + In 



+ ... 



(3 /2) 

This asymptotic relation will be of special value to our determination of i^g since it con- 
tains information on the vacuum matrix element (Os)^ (c/Eq. (|D). The large-s behavior of 
the spectral functions can be inferred from the logarithmic term in Eq. (|Tl|) via continuation 
to the real g^-axis. 



(11) 



{PV,3- PA,3){s) 



+ 



(12) 



Together, the spectral relations of Eq. ([ToD and Eq. ([12D imply 



Q%Uv,3-nA,3W 




+ r ds S^ [pv,3 - PA,3] (S) + OiQ-^) 

Jo 



(13) 



A. Correlators in d-Dimensions 

Consider the definition of Ily^s — 3 as expressed in d-dimensions, 

pt'^ J d'x e^'^-{Q\T{V,^{x)VaQ) - A|^(x)A;^(0)) |0) 

= {q^q'' - q'gn {^v,3 - H^.a) - q'q'TlZiQ') ■ (W) 

The energy scale p^.r. ('d.r.' denotes dimensional regularization) has been introduced to 
maintain the proper dimensions away from d = 4. It is straightforward to invert Eq. (p!4[) 
and we find 



3 



(0|T(V3''(x)\/„,3(0)-^;(x)A„,3(0))|0> 

Up to this point the procedure is well defined, as all quantities are finite-valued. 

To obtain a relation for {Oi)^, we need to evaluate Eq. (^) in the limit of x — 0. 
However, the asymptotic condition of Eq. ([Tl|) implies that unless the integral on the right- 
hand-side of Eq. (|15|) is regularized, it will diverge as x — 0. There are a number of ways 
to perform the regularization, and we shall consider two particularly useful approaches — 
first a momentum space cutoff directly below and then MS renormalization in Sect. IV. We 
shall distinguish vacuum matrix elements in the two schemes by means of the superscripts 
'(c.o.)' for momentum-cutoff and '(MS*)' for modified minimal subtraction. 



B. Two Sum Rules in Momentum-space Cutoff Renormalization 

Let us remove the divergence which occurs for = 4 in Eq. ( |15|) by cutting off the 
Q^-integral at = fi"^, where fi is the renormalization scale and for convenience we set 
fJ'd.r. = /U- It is valid to take = 4 in this case since the integral is finite. We find 

{Oi^r^ = ^ if dQ' Q\Uy,, - UA,s)iQ') ■ (16) 

Using Eq. (|T0|) to express this relation in terms of spectral functions, we arrive immediately 
at the following sum rule, 

^m^r^ = h^ f ds In [pv,3 - pas] is) . (17) 

It is equally straightforward to derive a sum rule for (a.Cg)},""-^- We first set 

= in 

Eq. (|llD to obtain 

(«.c)8)!r-°-) = ^(nv,3-n^,3)(p') . (18) 

Zn 

Because the variable is constrained in Eq. (0) to lie in the range where pQCD makes 
sense, the same must be true for the scale /i. Then by combining Eq. (|18D with Eq. (p!0D, we 
obtain the sum rule 

27r(a,08)f =h = r ds s^-^ [pv,3 - Pa.s] (s) . (19) 

Despite their apparent similarity, it is important to understand that there is a basic 
difference between the sum rules for {Oi)^^'°''' and {asOs)^^'°'^ ■ The former is obtained rather 
directly by taking the x — limit of Eq. ([15|) and using a cutoff in momentum to regularize 
the procedure. However, the latter rests upon assuming the dominance of the leading 
term in the OPE of Eq. (|Tl|) . This assumption becomes increasingly questionable as /x is 
lowered to energies just above the resonance region. It leads to an uncertainty in the value 
of {asOs)^^■°^^ which is not present in {Ol)^^■°^\ We postpone discussion to Sect. V regarding 
numerical evaluation of the integrals /i. Is appearing in Eqs. (|1^,(|19|). 
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III. KAON-TO-PION MATRIX ELEMENTS OF A LEFT-RIGHT OPERATOR 



A distinct but equivalent path to learn about {Oi)fj_ and (Cs)/^ is to perform a chiral 
analysis of the kaon-to-pion matrix elements themselves. However, the usual (V-A) x (V-A) 
weak hamiltonian 7iA5=i would be of no help in the chiral limit since its K-to-pi matrix 
elements vanish there. Instead we introduce a (V-A)x(V+A) nonleptonic operator (9a5=i 
defined as 



Cas=i 



9l 



-T 
2 

1 

-T 
2 



dix)ril + iXx) n(0)7"(l - 75)s(0) 



(20) 



where V^^, is the VT-boson propagator and V^, (a = 1, ... 8) are the flavor-octet vector, 
axialvector currents. Operators similar to Oa5=i have received some previous attention in 
the literature. The LR chiral structure of Oas=i ensures the survival of the K-to-pi 
matrix element Ai{p) = {7i~{p)\Oas=i\K~{p)) in the p —>■ limit, where we obtain 



M 



lim M (p) 

p=0 



9l 
16F2 



d^x V{x, M^) (0|T {V^^{x)V^M - A^(x)A^,3(0)) |0) . (21) 



A. Leading-log Analysis of QCD Corrections 

In the following, we perform a leading-log calculation of QCD corrections to the chiral 
matrix element M.. This leads naturally to renormalization group equations (RGE) for the 
quantities (Ci)^ and (ttsCs)^. 

Since the VT-boson propagator in Eq. ( pTj) acts as a cutoff for contributions with \x\ > 
we consider the leading term of the following operator product expansion, 

V,'{x)V,^{^) - A^,{x)Al{Q) = Vrm,^{0) - A^Am + 0{x) . (22) 
Evaluation of the spacetime integral in Eq. (^) is straightforward, 

Jd'xV,^{x,M'^) = ^ , (23) 

so that the matrix element specified at energy scale Mw becomes 

Mc.^^^{0,)m^ . (24) 

In order to express this vacuum matrix element at some lower energy fi, we must take 
QCD radiative corrections into account. The effect of these will be to mix Oi with Oq. The 
result of mixing at one-loop order is 
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47r 



3/2 
16/3 7 



0i' 

Os 



(25) 



where = qI/^t^ is the QCD fine structure constant and /i < My/. Using standard 
techniques we use the renormahzation group (RG) to provide a summation of the leading- 
log dependence over the range from Mw down to /x, 



M 



2V2F2 



[cM{0^)^ + Csi^^){Os),] 



(26) 



where 



Cl(/^) = g 



C8(/i) = g 



as(/i) 
,«s(M, 



8/9 



-1/9 



+ 



Q5(/^) 

,«s(M, 



8/9 



as(/i) 



a,(M, 



w) 



(27) 



with 



tts(/^) 



1 + 9^ In 
47r 



as(M, 



(28) 



An expansion of Eq. (^) through second order in as{fJ.) gives 

G p 



M 



+ 



3, /M|- 



327r2 



Stt 



-(a^0i),-(a,2O8), 



(29) 



Let us gain some feeling for the numbers involved. The minimum value of renormahzation 
scale considered in this paper is /iq = 2 GeV. Adopting this scale and taking as{Mw) = 0.119 
and as(/io) = 0.334 we find for the RG coefficients in Eq. (|26|) , 



M 



Gr 



2V2F^ 



1.071(Oi)^„ + 0.268(08) 



Mo 



(30) 



whereas the coefficients in the perturbative expression of Eq. ( [29| ) become 
Gp 



M 



{Oi),, + 0.88(«,O8)^„-0.52((«fO8)^„--(afOi)^„ 



(31) 



Finally, the condition of scale independence for the matrix element Ai, 



d 



/i 



(32) 



leads directly to the renormahzation group equations 



6 



d 



in 
1 



16 



(33) 
(34) 



To summarize — the above operator-product analysis involves computing radiative cor- 
rections perturbatively to one- loop order in QCD (c/Eq. (|25D ) and retaining only the de- 
pendence on leading logarithms in the evolution from scale M^y down to scale fi. The value 
of /i cannot be taken too small, otherwise the perturbative framework breaks down. 



B. Verification of the Operator Product Expansion 



Despite the explicit difference between the procedures of Sect. II and that carried out 
directly above, the two are equivalent. In particular, we can show that the {Oi)^ sum rule of 
Eq. ([l7|) and the (Os)^ sum rule of Eq. ( [l9|) reproduce the OPE to the leading log level. This 
verifies both the derivations that we provided and gives a direct insight into the workings of 
the OPE. 

Consider a partition of Ai characterized by the scale /i, 

M = M<{^i)+MM , (35) 

where A1<(/i) and A^>(/i) are dependent respectively on contributions with Q < ji and 
Q > H- Also, in addition to maintaining the requirement that /i lie in the pQCD domain, 
we further constrain it to obey /i <^ M]y. We then obtain 



MM 



32V27r2F2 



+ M2 



w 



dQ' ny,3(Q') - IIaAQ') + Oifx'/M, 



W) 



and 



X>(^) 



w 



dQ^ 



(36) 



(37) 



Upon inserting the large-Q form of Eq. ([IT]) into Eq. (0), we obtain 

Gp 



MM 



V2F^ 
3 

~32^ 



Alnf4^ 



In^ 



Comparison of Eq. (|29|) with Eq. (pSD yields the relation, 

G p 



MM 



(38) 



(39) 



2v^F2 

We see that the operators that we originally defined independently of the weak interaction 
are in fact the ones that appear in the Operator Product Expansion. Of course, this is to 
be expected, but it provies an explicit pedagogical demonstration of the nature of the OPE. 
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C. Sum Rules and RG Relations 



To complete the chain of logic, we demonstrate consistency of the spectral function 
sum rules for and (Og)^^ " -* with the corresponding renormalization group relations 

obtained previously from the operator product expansion (c/Eqs. (p3D,(|34D). The RG equa- 
tion for {Oi)^^'°'^ is immediately recovered upon differentiating the sum rule of Eq. (|r^) and 
making use of Eq. (O) , 



_d_ 
dii' 



3 /i 



ds 



s + fi^ 
,(c.o.) 



[PV,3 - PA,3] (s) 



in 



(40) 



To derive the RG relation for {Os)^^'°'' we start with the sum rule of Eq. (JT 



_d_ 
djj? 



(CO.) 



/i2 d 



ds 



2ti dfi'^ Jo s + /i' 



[PV,3 - PA,3] (s) 



(41) 



The integral in the above is seen to be fi^ (ny 3(/i^) — 11^ 3 (/i^)). We replace it using the 
aymptotic expression of Eq. (|T3|) to obtain 



P 



_d_ 



(CO.) 



In 



A2 



+ dss^ [pv,3 - PA,3] {S) + C(/i-2) 

Jo 

from which the RG relation of Eq. (B3) follows directly. 



(42) 



IV. MS RENORMALIZATION 

The work of the preceding sections was based on a momentum-space cutoff renormal- 
ization scheme, which is useful in yielding sum rules directly related to experimental data. 
At the same time, however, it is distinct from the more standard MS prescription. In this 
Section, we demonstrate how to relate the two approaches. 



A. Short Distance Analysis 

Let us reconsider Eq. (|TB|). We can show how the cutoff renormalization is related 
to dimensional regularization by keeping the high-Q^ part of the integral in Eq. ([T5| ) and 
analyzing it in terms of an e-expansion. We divide the integral into integration ranges below 
and above /i^. For the part of integral with below /i^, we can let d — 4 and recover 
exactly the cutoff integrand of Eq. (|16D, 

mf-^ = + /J dQ^ Q\iiy, - u^,)m (43) 

The asymptotic tail can be analysed in d dimensions. This introduces scheme dependence 
depending on which method is used to define Dirac algebra away from four dimensions, e.g. 
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the naive dimensional regularization (NDR) and t'Hooft-Veltman (HV) schemes in which 75 
is respectively anticommuting and commuting. We find 



-{d - 1)Q'' ■ (ny,3 - IIa,3){Q') = 27ras{Os), Q''-' 



l+{ds + -]e 



(44) 



where e = 4 — d and dg = (^scheme is associated with the loop integration and scheme- 
dependence. The values of dg in the NDR and HV schemes are 



dg 



-5/6 (NDR) 
1/6 (HV) . 



(45) 



The > fi^ integral can then be performed with the result 



^-, + ln4. + ^ + 24 



(46) 



The MS prescription is a subcase of dimensional regularization in which the terms 2/ (4 — 
d) — 'y + In Air in Eq. (|46|) are removed in the renormalization procedure. This gives our 
desired relation in a given scheme, 



3a, /3 



in \2 



7^ + 24 {Os) 



(47) 



To derive an analogous relation between (Cs)^^ and (Cs)!^'" '' we employ the leading 
behavior of correlators and spectral functions in the MS renormalization prescription, which 
has been calculated using the NDR scheme, M 



and 



+ 



119 



12 



(48) 



{PV,3 - PA,3)(s) ~ 



+ ... . 



(49) 



Then by setting Q = /i in Eq. (^) and combining the result with Eq. (^91) we find 

-2 r/ llQa 



ds 



s + fl^ 



[pv,3 - PA,3] (s) = 27ra, 



1 + 



247r 



(50) 



Comparison with Eq. (^7]) then implies that the NDR matrix element is given to first order 
in as by 



119a, 
247r 



(Os) 



(CO.) 



a.c 



TC 



(51) 
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(a) (6) (c) (d) 

FIG. 1. Some QCD corrections: full theory (a)-(b), effective theory (c)-(d). 



B. MS Matching at One Loop 



In this section we perform the matching at one loop and verify the scheme independence 
of the result. The effective operator Oa5=15 is expressed in terms of local operators Oi and 



O 



(eff) 
AS=l 



G p 



[ci(^)C>i + C8(/i)08 



(52) 



Determination of the coefficients ci and cs proceeds in two steps: first calculate QCD radia- 
tive corrections in both the full and effective theories, and then 'match' the two calculations. 
We shall carry out this procedure at one-loop order of the QCD radiative corrections. Since 
ci = 1 + 0{a1) and cs = 0{as), this will yield a determination of cs. For definiteness, we 
consider the free scattering of zero momentum quarks and adopt a common quark mass m 
to serve as the infrared cutoff. 

In the full theory, evaluations of the one-loop radiative corrections like those displayed 
in Figs. |I|(a)-(b) are finite and yield 



M 



G p 



(gg|Ci|gg)tree + 



1] {qq\Os\qq)tree + Oial 



(53) 



The analogous calculation in the effective theory is divergent and must be regularized. 
We employ dimensional regularization which introduces the scheme dependence mentioned 
above. Our calculation of amplitudes like those in Figs. [|(c)-(d) gives 



M 



G p 
2V2 



{qq\Oi\qq)tree + 



3a., 



Stt 

+ Cs{qq\Os\qq) tree 



1 + dse + dee) 



2 , . , 

7 + m 47r — m — ; 



{qq\Os\qq) 



tree 



(54) 



After the removal of the divergent term in MS renormalization, we compare the full theory 
and the effective theory to identify the coefficient function as 



In 



w 



--2d, 
2 



(55) 



When C^^gLi is applied to our problem of vacuum matrix elements, we have the amplitude 

Gp 



M 



272^2 



Cl(^)(Ol)r' + C8(/i)(08) 



MS 



(56) 
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Given our previous identification of tlie sclieme dependent operator Oi in Eq. (|47|), it can 
be seen tliat tlie sclieme dependence cancels between that of the matrix element and of the 
coefficient in the operator product expansion. 



V. NUMERICAL ANALYSIS AND UNCERTAINTIES 

We base our numerical determination of (Og)^^ "'-* and {Oi)^^'"'^ respectively on the sum 
rules in Eq. (|T9|) and Eq. (p!7[). This involves calculation of the integrals /§ and Ji, which 
contain the combination of spectral functions pv,3 — Pa,3, 

POO 

Ii= dsKi{s,p^) [py.3 - Pa.s] (s) (i = l,8) , (57) 
Jo 



with 



K, = s'^—^, K,^s'\r.r—^\ . (58) 

S + /i^ \ s ) 

As such, Jg and /i belong to the family of spectral integrals which include the DM0 sum 
rule H] , the two Weinberg sum rules [0] and the sum rule for the pion electromagnetic mass 



splitting The kernels occurring in these 'classical' sum rules are 

^DMO = , ^Wl = 1 , ^W2 = S , Xem = S In S . (59) 

This happenstance is most fortunate as the integrals defined by the kernels in Eq. ( p9D form 
a powerful set of constraints for any evaluation of /§ and Ji. Using an updated form of our 
earlier study ITT]] of chiral sum rules, we find for renormalization scale fi = 2 GeV the values 



Jg = -(0.30 ± 0.04) ■ 10-2 , Ji = -(0.42 ±0.06) ■ 10-2 , (60) 
At the higher renormalization scale, /i = 4 GeV, we obtain 

Jg = -(0.43 ± 0.06) • 10-2 , Ji = -(0.97 ±0.12) ■ 10-2 . (61) 



A. Uncertainties from Data Analysis 

The error bars quoted above correspond to our estimate of the uncertainty in the sum 
rules due to imprecision in our knowledge of the spectral functions. Before proceeding let us 
describe how these were arrived at, and assess other sources of uncertainty. The data at lower 
values of s are extremely well known, and they introduce very little uncertainty compared 
to other sources which we are concerned with. The high energy tail of py — Pa is small above 
s = 5 GeV2. In the p = 2 GeV integrals, there remains essentially no sensitivity to this high 
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energy tail once the constraints are imposed^. It is in the matching of these two regions 
that one encounters the greatest uncertainties. Fortunately, the four integral constraints 
described above are very stringent and allow us to limit the uncertainties in this region. 
We have used several methods to construct spectral functions which match the data within 
error bars and yet satisfy our sum rule constraints. These give variations in our integrals 
of under 6%. In addition, we have considered the situation where the asymptotic form of 
the spectral function is reached only on the average, with a damped oscillating term that 
provides deviations from the average. Since our sum rules are equivalent to transforming 
back to euclidiean Q^, these oscillations give exponentially suppressed effects at large /x^ 
once integrated. Again the constraints are very powerful in further limiting this effect, and 
our studies lead us to increase the uncertainty in the fit to 10% to account for this form of 
variation. 

We also must account for the fact that the data and the input into the constraints are 
measured in a world where is not zero, yet we are interested in the result in the chiral 
limit. This introduces corrections of order m^/m^ which is of order 3%. In fact, since we 
know some of the physics involved in passing to the chiral limit, we could attempt to make 
a realistic correction for the extrapolation to the chiral limit. However, since this would 
appear to introduce some model dependence into our procedure, we prefer to simply include 
the uncertainty as an error bar. In practice, the effect which has the most sensitivity for our 
results is the constraint of the pion electromagnetic mass difference, since the kernel i^'em 
bears the greatest resemblance to Ki and K^. Work on the pion and kaon electromagnetic 
mass differences indicates that the quark mass corrections are somewhat larger than average. 
Therefore, to be conservative we triple the canonical error estimate, leading us to quote a 
9% uncertainty for the extrapolation to the chiral limit. We have added this in quadrature 
to the statistical uncertainty to arrive at the error bars cited above. 

B. Uncertainty from the Operator Product Expansion 

Finally, we need to address the fact that it has become common to cite matrix elements 
at a scale = 2 GeV, which is a rather low scale for perturbative QCD to be fully in the 
asymptotic region. In fact, our method can be used for any /i, and we can check if the 
asymptotic QCD behavior is obtained. For example, the renormalization group equations 
relate the //-dependence to the magnitudes of the operator matrix elements. Although one 
of the relations (Eqs. (^),(^)) is automatically satisfied, we explicitly showed above that 
the second holds only if /i is large enough, i.e. that it is well into the region where the 
asymptotic tail of the spectral functions becomes applicable. It is easy to see from the data 
alone that this is not the case at // = 2 GeV. Another way to state the same result is that 



^At higher values of /i there occurs more sensitivity to the asymptotic tail, and it is the tail that 
describes the logarithmic running of the Og matrix element. While we have a good handle on the 
size of the 1/s^ component of the tail, we know little about the \/s^ component. However as long 
as the portion is not much larger than the behavior for s > 5 GeV^, its effect is within 
our quoted error bars. 
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there remain power corrections in the sum rule, although the renormalization group states 
that the running with fi should be only logarithmic. We believe (because of the generality 
of our framework) that this issue must also be present in the lattice results, and we urge the 
evaluation of lattice matrix elements at larger values of /i. 

We do see such non-asymptotic behavior in our results. Actually, for Oi our method of 
cutting off the high frequency modes of the current-current product is in accord with Wilson's 
original idea of the definition of a scale-dependent matrix element. Therefore, our sum rule 
for can be treated as a definition of this amplitude at any scale fi, even if that scale 

is not yet asymptotic. For Og, however, there is some uncertainty as to an ideal definition in 
the non-asymptotic region. For example, the RG relation of Eqs. (p3D,(^0D requires that we 
use exactly our definition, yet this only is foolproof if the RGEs are fully valid. Equivalently, 
if this matrix element is defined via the coefficient of in the vacuum polarization, there 
can be order corrections remaining if one works in the non-asymptotic region. We see 
evidence of such power corrections. Moreover, attempting to discard the effects leads to 
a larger matrix element. At /i = 4 GeV, the corrections are rather modest, in line with other 
uncertainties that we have described. However, at /i = 2 GeV, these non-asymptotic effects 
represent a significant intrinsic uncertainty. We have taken these into account by combining 
two evaluations, one obtained by evaluating the sum rule at fi = 4 GeV and using the RGE 
to transform down to n = 2 GeV, and the other by direct evaluation of the sum rule at the 
lower scale. We average these two and assign the difference as an independent error bar. 
The error bar is chosen such that a one-sigma variation reproduces the full range between 
the two methods of evaluation. We do this for both matrix elements. The quoted error bar 
at /i = 4 GeV is scaled down from the fi = 2 GeV values by a factor of four, as appropriate 
for quadratic power corrections. 



C. Conversion to MS Renormalization 

We now transform to the MS matrix elements. The results of our direct evaluation at 
fi = 2 GeV leads to the matrix elements 

((^8)& = -(0.67 ± 0.09) ■ 10-3 GeV*^ , 

{<^i)?^Iy = -(0-70 ± 0.10) ■ 10-^ GeV*^ , (62) 

where we have taken as{2 GeV) ~ 0.334. When we evaluate the integrals at ;U = 4 GeV and 
use the RGE to rescale back to /i = 2 GeV, we instead obtain 

{Os)f^ly = -(1.29 ± 0.15) ■ 10-=^ GeV*^ , 

= -(1-02 ± 0.10) • 10-^ GeV^ , (63) 

which is a measure of the potential non-asymptotic corrections found at low values of fi. As 
explained above, this leads us to quote our result as 

Wf^lw = -(0-98 ± 0.13 ± 0.23) ■ lO^^ GeV*' , 

(Ci)2'1>v = -(0-86 ±0.10 ±0.16) • 10-^ GeV'^ , (64) 
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The first error bar corresponds to the uncertainty in the evaluation of the sum rule whereas 
the second is the potential non-asymptotic uncertainty defined above. Note that the two 
matrix elements differ by an order of magnitude. The related K'^ vrvr matrix elements are 
then 

((7r7r),=2|Q?^'^|i^°)2 GeV = (0.43 ± 0.05 ± 0.08) GeV^ , 
((7r7r),=2|Q?^'^|i^°)2 GeV = (2.58 ± 0.37 ± 0.47) GeV=^ . (65) 

In the NDR scheme with fi = 2 GeV this translates into the following B-factor determina- 
tions, 

5f/'^[NDR,/i = 2 GeV] f^^^-^) = 0.55 ± 0.07 ± 0.10 , 

Knis + m J 

Bi^^'^\NDR, = 2 GeV] ( ^'^ ^^^ ) = 1.11 ± 0.16 ± 0.23 . (66) 

Note that the combination of B-factor and quark masses is 'physical', appearing in the 
formula for e'/e. The comparison of these results with some lattice evaluations is hampered 
by the fact that our evaluation is of the full matrix elements, while most lattice calculations 
are of the B- factors directly [|12[. If large values of quark masses are used, our results are 
larger that other estimates, yet for the currently favored smaller quark masses the results 
are not inconsistent. There is one recent lattice evaluation which provides absolute matrix 
elements which we can compare to. The Rome group fl^ quotes the matrix element in 
in the the quenched approximation using K ^ it matrix elements plus the chiral relation 
between K ^ n and K ^ titi. When the meson masses are taken as the kaon mass they 
find, in the NDR scheme at = 2 GeV, 

{{7in)i=2\Q?'^\K')2 GeV = (0.22 ± 0.04) GeV^ , 

{{7i7i)j=2\Qf^'^\K')2 GeV = (1-02 ± 0.10) GeV^ . (67) 

The quoted error does not include estimates of the effect of quenching nor the extrapolation 
to the continuum limit. Their results seem to be systematically smaller than ours. 
Finally at /x = 4 GeV we have the vacuum matrix elements, 

(C8)ffev = -(1-63 ± 0.20 ± 0.06) ■ 10-3 GeV^ , 

(C'i)fSv = -(1-71 ±0.20 ±0.04) ■ 10-^ GeV^ . (68) 
The corresponding K —>■ tttc matrix elements are 

((7r7r)/=2|Q?^^^|^°)4 GeV = (0.85 ±0.10 ±0.02) GeV^ (69) 
((7r7r)/=2|Q8^^'^|i^°)4 GeV = (4.34 ± 0.56 ± 0.15) GeV^ 
and for the B-factors we obtain, 

sJ'/')[NDR,/i = 4 GeV] (^^-^^ = 1.10 ± 0.13 ± 0.03 , 

V ms + m J 

B^s^^\nDR,h = 4 GeV] ( ^'^ ^^^ ) = 1.87 ± 0.25 ± 0.07 . (70) 
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That B^^"^^ has a large variation with is expected from the RGE of Eq. (|34D , given our 



previous resuh that the vacuum matrix element of is much larger than that of 0\. 

VI. CONCLUDING COMMENTS 

The method that we have described has the virtue of being a fully rigorous framework. 
Moreover, the input data is largely taken from experiment, and hence represents an evalua- 
tion that is model independent. Besides the direct comparison with the results with lattice 
calculations, there may also be other lessons in this calculation. Since in our method the 
matix elements are evaluated by constructing the Euclidean vacuum polarization function, 
lattice calculations may also be able to directly follow many of the steps in our procedure, 
and thereby test their methods in more detail. By explicitly studying the product of currents 
at non-zero values of the spatial separation, the matrix elements can be evaluated without 
some of the operator mixing problems that occur on the lattice when using local operators. 
Moreover, by studying vacuum matrix elements as well as hadronic matrix elements, the 
chiral relations can be checked on the lattice. Finally, we recall the lesson, discussed above, 
that power-law corrections still appear to exist at /i = 2 GeV, especially in the Og matrix 
element. This raises the concern that when one is working at such a low value of /i, there 
may be significant corrections even in lattice evaluations. Certainly, use of ;U < 1 GeV, as 
occurs in many model dependent evaluations, appears extremely dubious. 

The values displayed above are based on working in the chiral limit of massless quarks. 
One must, however, add to these the chiral corrections. Work has begun on this important 
problem. Nevertheless, it is interesting to look at the phenomenological consequences of 
the results reported in this paper.. While we cannot give a full evaluation of e'/e because 
we have not evaluated the contribution of -Bg, we can give the contribution arising from the 
electroweak penguin, 

= (-12 ±3) ■ 10"^ . (71) 

The effect of is expected to be positive, and needs to be almost three times larger than 
that of i?8 if the Standard Model is to account for the experimental value. 
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